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^ We study the potential of the LHC and future low-energy experiments to precisely 

measure the underlying model parameters of a new Z' boson. We emphasize the com- 
plimentary information obtained from both on- and off-peak LHC dilepton data, from 
the future Q-weak measurement of the weak charge of the proton, and from a proposed 
measurement of parity violation in low-energy M0ller scattering. We demonstrate the 
4^ importance of off-peak LHC data and Q-weak for removing sign degeneracies between 

^ Z' couplings that occur if only on-peak LHC data is studied. A future precision mea- 

surement of low-energy M0ller scattering can resolve a scaling degeneracy between 
^ quark and lepton couplings that remains after analyzing LHC dilepton data, permit- 

'— ' ting an extraction of the individual Z' couplings rather than combinations of them. 

^ We study how precisely Z' properties can be extracted for LHC integrated luminosities 

^ ranging from a few inverse fcmtobarns to super-LHC values of an inverse attobarn. For 

the several example cases studied with Mz' = 1.5 TeV, we find that coupling combi- 
nations can be determined with relative uncertainties reaching ±30% with 30 fb~^ of 
integrated luminosity, while ±50% is possible with 10 fb~^. With SLHC luminosities 
\C) of 1 ab~^, we find that products of quark and lepton couplings can be probed to ±10%. 
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1 Introduction 



Z' gauge bosons arise in most constructions of physics beyond the Standard Model (SM). 
They appear in grand unified theories such as S'O(IO) [1] and £'(6) [2], in Little Higgs 
models [3], and in theories with extra space-time dimensions [4]. They often appear as 
messengers which connect the SM to hidden sectors, such as in Hidden Valley models [5]. 
Due to their pervasiveness in models of new physics, significant effort has been devoted to 
searching for them experimentally. Z' states that decay to lepton pairs have a simple, clean 
experimental signature and can easily be searched for at high-energy colliders. Current direct 
search limits from the Tevatron require the Z' mass to be greater than about 1 TeV when 
its couplings to SM fermions are identical to those of the Z boson [6]. 

Since they are so prevalent in models of new physics, and it worthwhile to study what 
could be learned from a Z' discovery at the LHC This also serves as a useful LHC benchmark 
for how well new-physics Lagrangian parameters can be experimentally determined. Since 
the experimental signature is clean and the QCD uncertainties have been studied and found 
to be fairly small [7-9], it is likely that the couplings of a discovered Z' can be studied 
with reasonable accuracy to probe the high scale theory that gave rise to it. Many studies 
of how to discover, identify, or measure Z' properties and couplings to SM particles have 
been performed [9-16]. Recently, we performed a detailed analysis using a fully differential 
next-to-leading order QCD simulation to quantify how well all possible Z' couplings could 
be extracted from on-peak data at the LHC assuming realistic QCD, parton distribution 
function (PDF), and statistical errors [8]. Our study found that four combinations that 
probe both the parity-symmetric and parity-violating couplings of the Z' could be probed, 
with some combinations measurable with roughly 15% precision assuming 100 fb"^ of data. 

This previous study contains several limitations. More statistical power is possible if 
data from below the Z' peak is utilized. Also, additional differentiation between various 
Z' couplings is possible with off-peak data. On-peak, the only combinations of couplings 
that appears are (g x e)^, where q and e denote arbitrary quark and lepton couplings. Two 
degeneracies exist in this expression: the relative signs between various g x e combinations 
are not determined, and the scaling q — > yq^e e/y, with y an arbitrary constant, leaves 
the expression unchanged. The first can be removed by Drell-Yan production of the Z' 
off peak at the LHC, where the cross section contains a linear dependence on g x e. The 
removal of the second scaling degeneracy requires observation of the Z' in a channel other 
than Drell-Yan at the LHC, or in a different experiment. 

In this manuscript we extend the previous study to address several of these issues for the 
example case of a Z' with mass Mz' = 1.5 TeV. We perform a detailed simulation of the 
Z' signal and SM Drell-Yan background to study the potential of the LHC for a precision 
extraction of the model parameters. We include errors arising from statistics, uncertainties 
in PDFs, and missing higher-order QCD and electroweak corrections. We study how well 
couplings can be measured assuming different energy and luminosity scenarios at the LHC, 
including collisions at both y/s — 10 and 14 TeV, and integrated luminosities ranging from 
3 fb"^ to super-LHC amounts of 1 ab ^. li a Z' is found at the LHC, additional information 
on its parameters can be obtained from various low-energy measurements. Much of this 



1 



information is complementary, as many low-energy observables depend on different coupling 
combinations than LHC measurements. We study several planned future experiments with 
the potential for observing and studying a deviation from their SM expectations, and that 
may help remove the degeneracies described above. We focus on the Q-weak experiment [17], 
which plans to measure with high precision the weak charge of the proton in electron-proton 
scattering, and a proposed Jefferson Laboratory measurement of the effective weak-mixing 
angle through parity violation in low-energy M0ller scattering [18]. This second experiment 
is of special interest. Since it depends only upon the leptonic Z' couplings, it can potentially 
break the scaling degeneracy between quark and lepton couplings to the Z' . Our primary 
findings are summarized below. 

• Discrimination between the studied test models at the 68% and 90% confidence level 
with as little as 3 fb~^ of LHC data at the center-of-mass energy a/s = 10 TeV. With 
10 fb~^ of integrated luminosity at this energy, measurements of the Cg parameters 
defined in Refs. [8, 15] with relative uncertainties of ±50% are possible, while ±30% 
errors become possible with 30 fb~^. 

• With high luminosities of 300 fb~^ at ^J~s = 14 TeV, the relative signs between the 
various q x e combinations are determined, and measurements of these combinations 
approaching 10% precision are possible. 

• An interesting feature we derive from our analysis is an independent extraction of the 
width Tz' to a few GeV precision obtained without fitting the Z' Breit-Wigner peak. 
Since the on-peak cross section scales as ~ 1/Tz' while the off-peak result contains no 
width dependence, comparison of data from different invariant mass bins allows a good 
determination. 

• A precision measurement of the weak charge of the proton with Q-weak helps remove 
sign degeneracies between the q x e that remain after LHC running. 

• With the inclusion of a future measurement of the effective weak-mixing angle in M0ller 
scattering, the individual q and e charges can be separately determined. We demon- 
strate for an example model how a combination of LHC Drell-Yan data and M0ller 
scattering make this possible. LHC measurements restrict the Z' model parameters to 
lie in a wedge in the versus coupling plane, while M0ller scattering restrict the 
couplings to a hyperbolic region. A combination of the two experiments localizes the 
couplings to a point in this plane, within errors. 

Our paper is organized as follows. In Section [2] we review the parameters which describe 
Z' production at the LHC and introduce the various test models we use to illustrate our 
analysis procedure. We describe the LHC observables which we use to analyze Z' properties 
in Section |3| and discuss the expected low-energy measurements of the proton weak charge 
and the effective weak-mixing angle in Section |4j We present our analysis of Z' couplings 
at the LHC and in the Q-weak experiment, which are both sensitive to the q x e coupling 
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combinations, in Section |5j We study both low and high integrated luminosities in this 
Section, detailing the expected uncertainties in the parameter extractions for a variety of 
different running scenarios. A summary showing how well all g x e couplings can be measured 



as a function of integrated luminosity is given in Section 5.6 In Section ^ we consider the 



constraints imposed by a proposed M0ller scattering experiment, and demonstrate how a 
precision measurement at J-Lab allows the determination of the individual q and e charges. 
Finally, we conclude in Section [7) 



2 Theoretical Framework 

In this section we motivate a parametrization of Z' models which can be probed by future 
experiments, while remaining as model-neutral as possible. We also describe examples of 
common models which fit into our parametrization, and which we use to test our procedure 
for extracting these parameters. 



2.1 Z' parameters 

At the LHC, the smoking-gun signal for a is a new resonance in dileptons. Other particles 
that decay into such final states, such as Kaluza-Klein gravitons, can be distinguished from 
a Z' by study of the lepton angular distributions [19,20]. The production of this resonance 
depends on the following parameters: 

• the location of the resonance, i.e. the mass of the state, M^/; 

• the width of the resonance, T z''-, 

• the couplings of the quarks which produce the Z', u'^j^,u^j^,d^]^,d^^ with the index i 
denoting the generatioiiQ 

• the couplings of the charged leptons in the observed decay mode, e^, e^, and the 
neutrino coupling v\^. 

While the Z' may couple to and decay into any number of other particles, from W bosons 
to new exotic fermions, these parameters only enter into the Drell-Yan channel we consider 
through the width. 

One cannot hope to determine all of these parameters in a completely model-independent 
way. To proceed, we must make some assumptions on the parameter space. Fortunately, sev- 
eral restrictions are very well motivated by experimental data. We first make the assumption 
that the Z' couplings are generation independent. If they were not, a Z' light enough to be 
discovered at the LHC would generate large flavor changing neutral currents in contradiction 
with experiment. Z' bosons with coupling strengths approximately the same as the SM Z and 
(9(1) coupling differences between the first two generations must have masses greater than 

^We have absorbed any overall coupling into these parameters to simplify the discussion and analysis. 
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roughly 10^ TeV due to constraints on the generated AF — 2 four-fermion operators [21]. 
Constraints on Z' bosons with different couphngs to third-generation fermions are weaker; 
deviations from our parameterization for top quarks, bottom quarks, or tau leptons can be 
searched for explicitly by observation of Z' decay into these final states [22-24]. Second, we 
assign states in a given SU{2)l doublet the same coupling to the Z' . If the generator of 
the new gauge group to which the Z' belongs does not commute with SU{2)l, the Z' would 
gcnerically couple to the electrowcak symmetry-breaking sector, inducing Z — Z' mixing. 
Specific models which would protect against this mixing, such as Kaluza-Klein Z's, should 
be distinguishable from the class of models considered here [14]. LEP Z-pole measurements 
restrict the mixing angle to be smaller than a few thousandths of a radian [25,26]. We take 
this magnitude of mixing to be negligible. We note that arranging the Z-Z' mixing angle to 
be this small requires, for the models considered, either fine-tuning the charges and vevs in 
the Higgs sector which gives mass to the Z and Z' bosons, or setting these parameters smaller 
than their expected values by roughly an order of magnitude [15]. There are many models 
in the literature which avoid these issues, but the Eq models we use as examples remain a 
useful benchmark, particularly for comparing with other studies. These restrictions leave 
the following seven Z' parameters to be determined: 

• the mass Mz' and the width Vz'] 

• cl and cr, the couphngs to the lepton doublet and right-handed electron; 

• g^, the quark doublet coupling; 

• ur and dij, the right-handed quark couplings. 



To illustrate how well our analysis strategy determines Z' properties, we test it on several 
example models. Several well-known grand unified theories (GUTs) fit into the framework 
described above. We take three such examples to illustrate how well their parameters could 
be determined if they are found. Two of these arise from the exceptional group Eq^ and one 
from 5*0(10). Below we briefly describe these models and list their couplings. For a more 
comprehensive introduction to these specific examples, we refer the reader to several excellent 
reviews [1,2]. Our procedure is by no means limited by these choices; we choose these 
common examples simply to demonstrate the efficacy of our model-independent approach. 
This method can be used for any Z' effective theory satisfying the coupling assumptions 
outlined above. 

• Eq Models : Eq models are described by the breaking chain 



2.2 



Test models 



Eq 



50(10) X U{1)^ SU{5) X U{1)^ X U{1)^ ^SM X U{1)^ 



(1) 



where 



Z' = Z'^ cos p + Z'^ sin p 



(2) 
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is the lightest new boson arising from this breaking. In this paper we take the x model 
(/? = 0), and the ifj model (/5 = 7r/2), as representatives, though the whole family fits 
into our framework. 

• Left-right models: We also consider a left-right model coming from the symmetry group 
SU{2)fi X SU{2)l X U{1)b-l- Left-right models can arise from the following breaking 
of SO (10): 

SO{10) ^ SU{3) X SU{2)l X SU{2)n x U{1)b-l. (3) 
The Z' in left-right models couples to the current 

Jm = c^lrJTr - ll1a^RJT-L. (4) 

with aLR= y^icwdR/ ^w9l) ~^ 1' ^^'^ gi = e/ cos 6w- We examine the symmetric case 
Ql = dR, where ulr — 1-59 if one takes the on-shell value of sin^ 9w- 

The fermionic couplings of these models are summarized in Table [T] For the standard 
definition of the Eq models, we take the overall coupling to have its GUT-scale relation to the 
EM coupling down to the Z' scale. An overall e/ cos 9w has been factored out. We will later 
also consider a version of the x model with a slightly larger overall coupling to demonstrate 
what is gained from potential measurements in low-energy experiments. The overall factor 
of e/ cos6'vi/ — 0.36 in this model has been replaced by 1/2. We denote this model x*- 
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-1 
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12 
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2v^ 


-v^ 
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-3 
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-v^ 
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3 
2v^ 


12 


1 

2oLij 


dR 


1 

2v^ 


12 


1 "Lfl 
2aLR 2 



Table 1: Fermion couplings to the Z' for the considered models. An overall e/ cos 9y/ has 
been factored out. 



3 LHC Analysis Framework 

The primary means by which we will probe Z' properties will be through examining the 
dilepton mode at the LHC. We describe here the structure of the Z' cross section, present 
the observables we utilize for various choices of integrated luminosity, and discuss the details 
of our simulation procedure. 
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3.1 LHC measurements 



We would like to determine seven parameters: M^', Tz', Ql, ur, da, cl, and e^. The location 
of the resonance peak in the invariant mass of the dileptons should determine M'^ extremely 
well with a sufficient number of events; we take the error in this measurement to be negligible. 
For this analysis, we take as an example Mz' = 1.5 TeV. For most models, much lighter Z's 
are excluded by experiment [6,25]. For significantly heavier Z's, testing shows that there 
is not much improvement over our previous on-peak analysis, [8]. In principle, the width 
can be determined from the resonance shape if the experimental resolution in invariant mass 
does not dominate the natural width of the Z'. With enough statistics, widths smaller than 
the experimental resolution could possibly be probed by fitting the observed spectrum to a 
convolution of Z' resonance shape with the experimental Gaussian resolution. How well this 
can be done requires a detailed detector-level study and is beyond the scope of this analysis. 
We make no assumption on how well the width can be determined, beyond the requirement 
that one be able to define an "on-peak" bin containing most of the resonance, to be explained 
below. For an experimental study of both the Z' mass and width measurement at the LHC, 
see Ref. [27]. 

One can reconstruct the following kinematic variables from the dilepton measurements 
at the LHC. Taking the four-momenta of the electron and positron to be ki and k2, we have, 
with kz' — ki + k2: 

• the invariant mass of the pair, — kzi] 

• the reconstructed rapidity of the Z', Y — ^ log ^' ; 

• the scattering angle of the electron, measured in the CoUins-Soper frame [28]. The 
sign of cos^* is chosen to be positive if k\ and A';|, have the same sign, i.e., the electron 
is scattered in the same direction as the Z' is boosted down the beam pipe. 

The differential cross section of pp — > e^e" can then be written in the following form: 
dYdcZ*dM^ ^ ^ ^^'^^^^ + ^^^^'^ + + 4iMz',rz'){ql - ql){el - el) 

q=u,d 

+bl{Mz')qReL + bl{Mz')qReR] 

+h{Mz')qLeL + b4{Mz')qLeR + c. (5) 

The coefficients a, b, and c contain all kinematic and PDF dependence. The c term is the 
contribution from the SM 7 and Z. To predict a measurement, one only has to integrate 
them over the appropriate region, and then input the dependence on the Z' couplings. By 
choosing the regions to integrate over (measurement bins) appropriately, one can solve the 
system of equations and gain coupling information. We now motivate the choice of bins. 

• Forwards/Backwards: At the parton level, terms symmetric in the quark-lepton scat- 
tering angle go as {qRqji + q^qDi^^R^^R + e\e'^), where i, j e 7, Z, Z' . Terms antisym- 
metric in the scattering angle go as (g^g^ — ?L?i)(e^e]j — e^ej^). We could separate 
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these terms by binning events depending on whether they are forwards or backwards 
scattered; F + B is symmetric and F — B is antisymmetric. We do not have direct 
access to the partons, however, and the LHC is a pp colhder, so there is no preferred 
direction for the quarks. We define "forwards" to be F = f^d cos 9* and "backwards" 
to be S = J^^ dcos 9*. The advantage of defining the direction in terms of the Z' boost 
is that this direction correlates with the quark direction, and therefore one still gets 
useful information from the asymmetry even though the LHC is a pp collider [29]. 

• Rapidity: The higher the rapidity Z', the more likely it is to have come from a valence 
u quark and sea u quark. This is due to the dominance of valence u at high x. In 
principle, to separate u from d-type quarks, one needs only two rapidity bins, where 
the dividing line is chosen such that the bins contain approximately the same number 
of events [11]. There is no penalty to choosing additional bins, as long as one does not 
bin so finely that detector resolution becomes an important issue. For the analysis with 
Mz' = 1.5 TeV, we choose the bins < |r| < 0.4, 0.4 < |r| < 0.8, 0.8 < |r| < 1.2, 
and 1.2 < |F| < Ymax, where Ymax = 1/2 log(s/M|,) is the largest rapidity available. 
We later study the effect of the coarser two-bin analysis, and find that only slight 
improvement is obtained by this finer binning. 

• Invariant Mass : Previous studies have examined the above observables and simulated 
extraction of coupling information at the LHC [8,11-13], or bounded couplings at the 
Tevatron [15], by integrating the invariant mass around the Z' resonance peak. Here, 
the b and c terms in Eq. ([s]) can be safely ignored, but any information they might 
provide is discarded also. Only squares of couplings are accessible on the resonance 
peak. We aim to gain sign information by probing regions where the Z' interference 
terms with the other gauge bosons have an effect, i.e., we want to study the b terms. 
We will do this by including invariant mass bins between the Z pole and the Z' pole. In 
the 1.5 TeV case, we take the bins 800 GeV < M < 1000 GeV, 1000 GeV < M < 1200 
GeV, 1200 GeV < M < 1400 GeV, and an on-peak bin. The on-peak bin should be 
chosen to contain most of the resonance; three or more widths around the peak should 
be sufficient. However, since we assume no knowledge of the width, we must fix a 
number beforehand. We have chosen the on-peak bin 1400 GeV < M < 1600 GeV, 
but the exact size does not affect the quality of the analysis, so long as it is consistent. 
For simplicity we refer to the bin as on-peak henceforth. The appropriate size to use for 
the on-peak bin can be determined with only a rough idea of the Z' resonance shape. 
The other bins should be chosen far enough away from the pole that the a terms do 
not dominate. For our considered models and mass, observables in each invariant mass 
bin receive contributions roughly equal in magnitude from the a, b, and c terms, with 
c being somewhat stronger in low mass bins and a being stronger in high mass bins. 

In total, we simulate 2 x 4 x 4 = 32 measurements for the main part of our LHC analysis. 
A summary of the bins used for Y, M, and cos 9* are shown in Table. [2j For illustrative 
purposes a selection of numerical results for both the SM and the test models using the 
factorization and renormalization scale choices fip = fi^ = Mz' can be found in Table [3j 
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Bin 1 


Bin 2 


Bin 3 


Bin 4 




[0,0.4] 


[0.4,0.8] 


[0.8,1.2] 


[1.2, Yfjiax 


M (TeV) 


[0.8,1.0] 


[1.0,1.2] 


[1.2,1.4] 


on-peak 


cos^* 


[-1,0] 


[0,1] 







Table 2: Summary of bin cndpoints used in our analysis for the studied variables: the rapidity 
Y, the invariant mass A/, and the leptonic scattering angle cos6*. Each region is displayed 
in the format [lower, upper]. Further details are given in the text. 





SM 


X 


LR 


i, 


X* 


(^on—peak 


0.783 


52.6 


59.7 


26.6 


94.8 


F 


< F < 0.4 


800 < M < 1000 


1.57 


1.68 


1.45 


1.59 


1.85 


B 


0.4 < y < 0.8 


1000 < M < 1200 


0.382 


0.524 


0.464 


0.378 


0.772 


F 


0.8 < Y < 1.2 


1200 < M < 1 100 


0.281 


0.395 


0.120 


0.383 


0.696 


F 


Y > 1.2 


on-peak 


0.075 


1.90 


3.92 


1.59 


3.39 


B 


Y > 1.2 


on-peak 


0.025 


2.95 


2.52 


1.52 


5.33 


F 


< F < 0.4 


on-peak 


0.147 


9.60 


11.34 


4.65 


17.29 


B 


< F < 0.4 


on-peak 


0.114 


10.65 


10.24 


4.61 


19.22 



Table 3: Selected expected cross sections for the considered models and the Standard Model, 
in fb. Parameters and selection cuts are discussed in the text. 

3.2 Simulation details 

To simulate observation of a Z' at the LHC we perform a fully differential next-to-leading- 
order (NLO) QCD calculation of all observables described above. We use CTEQ 6.6 PDFs [30]. 

The QCD corrections to the cross section are needed to obtain the proper normalization of 
the observables [8]. We impose the following basic acceptance cuts on the final state lepton 
transverse momenta and pseudorapidities: pip > 20 GeV and [r^'l < 2.5. Together, these con- 
straints result in acceptances of roughly 90% [8]. Previous studies have found that detector 
resolution effects and other measurement errors are unlikely to have a significant effect on 
the e+e" final state [13], and are neglected. A CMS simulation of Z' production found recon- 
struction efficiencies near 90% in the electron channel and no significant detector systematic 
errors [31]. In addition, electron energies can be measured to better than 1% accuracy, and 
invariant masses can therefore be reconstructed very well. Reconstruction efficiencies above 
90% and invariant mass measurements to sub-1% precision have also been found for high 
invariant-mass muon production at the LHC [32] . We include here only the statistics for the 
e"'"e~ production, without reconstruction efficiencies imposed. In light of the possibility for a 
factor of two increase in statistics from precision measurements in the muon-pair channel, we 
beheve our study conservatively estimates the LHC potential. Equivalently, one can rescale 
all required luminosities by a factor of 2 if including this channel, or use it to separately 
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extract electron couplings and muon couplings. 

To construct these measurements, we evaluate Eq. ([s]) in the following manner. The 
considered models have narrow widths. For computational efficiency, we calculate the a 
coefficients once using the x width of 17.9 GeV. We choose this width because it is centrally 
located among those considered. On the resonance peak, observables scale as l/F in the 
narrow-width approximation. Since we are not integrating over all invariant masses, we 
must improve this approximation for a finite bin size denoted B. When we scan over the 
width, we scale the a coefficients integrated over the on-peak bins by replacing 

/ rfM2a(r) (1 + 2^^) / dMMT^) (6) 

to integrate a(r) for other widths. Here we have included an 0{rz'/B) correction, and 
dropped terms of 0{B^/Ml). For M = 1.5 TeV, we have taken B = 200 GeV. The 
remaining integral of a(T^) is computed using the actual propagator, not with the narrow- 
width approximation. This is accurate to within tenths of a percent over the range of widths 
considered, and allows the a coefficients to be computed once rather than for all widths 
scanned over. If the finite bin-size corrections are completely dropped, errors of the size 
5 — 10% are introduced. We can then calculate the integrated coefficients a, b, c once the 
mass is known. We take c to be the SM Drell-Yan cross section. 



3.3 Early LHC running 

We will also consider subsets of the above measurements for use in discriminating models 
with early LHC data assuming = 10 TeV. Our restriction in the number of bins is moti- 
vated by the desire to avoid low-count statistical issues and the need to account for detector 
miscalibration, which would over-complicate our analysis. At low luminosities the reduc- 
tion in bins does not sharply reduce our sensitivity to differences in the underlying model 
parameters. We merge the bins as follows for 10 TeV running. 

• 30, 100 fb~^: The two highest off-peak invariant mass bins are merged (now 1000 < 
M < 1400), as are the two highest rapidity bins (now Y > 0.8), for 2 x 3 x 3 = 18 total 
bins. While the LHC will hopefully increase in energy prior to reaching this integrated 
luminosity, we consider these measurements at 10 TeV as a bridge between the early 
running and the late/SLHC running. 

• 10 fb~^ : All off-peak invariant mass bins are merged (800 < M < 1400), as are the 
two highest rapidity bins, for 2 x 3 x 2 = 12 total bins. 

• For the pairwise model comparison at a few femtobarns of integrated luminosity to 
be described in Section [5| we consider only measurements on the resonance peak. 
For simplicity we use four bins: 2 rapidity bins (|y| < 0.6, \Y\ > 0.6), and the for- 
wards/backwards measurement. 
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Below we orient the reader by providing 10 TeV on-peak cross sections for our test models. 
The parton luminosity for Mz' = 1.5 TeV is consistently reduced by approximately a factor 
of 2.5 from = 14 TeV. 



Model 


a (fb) 


X 


20.5 


LR 


23.8 




10.7 


X* 


36.9 



Table 4: On-peak cross sections for test models at ^/s = 10 TeV. 



4 Low-energy Measurements 

The upgrade of the Jefferson Lab polarized e~ beam will advance the precision frontier 
for low-energy parity violation. The Q-weak experiment [17] will probe the effective weak 
charge of the proton in ep scattering, and the proposed M0ller scattering experiment [18] can 
determine sin^ 9'^-^ through a precision measurement of the asymmetry in e~e~ scattering. 
Deviations from predicted SM values due to a Z' can provide further coupling information 
once the Z' mass is known from the LHC We describe here which Z' parameters to which 
each experiment is sensitive, and review the expected errors on the measured quantities 
arising from both experimental issues and SM theoretical uncertainties. 



4.1 Q-weak 

The Q-weak experiment will probe the effective weak charge of the proton, Q^^ ~ 1 — 
4 sin^ 9\Y, by observing the parity- violating asymmetry in ep scattering. Due to the smallness 
of 1 — 4sin^^vF! higher-order corrections are important. The predicted low-energy value in 
the SM (from running from Z-pole measurements) is [33] 

QP^ = 0.0713 ± 0.0008. (7) 

The anticipated experimental error after Q-weak running is 4.3% [17]. The experimentally 
observed value will differ due to the presence of a Z'; the shift is [34] 



^'^w = TTTT^C^^aUv + eAdv), (8) 



2V2 

where A and V denote the axial-vector and vector charges, respectively. We rewrite this in 
terms of the parameters defined in Section |2.1| to derive 



a/2 

^Qw = 2^2^^^ (3g^ + + dR)ieR - cl). (9) 
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Like LHC observables in the Drell-Yan channel, the Q-weak measurement is sensitive only 
to the coupling combination q x e. 



4.2 M0ller scattering 

The goal of the proposed J-Lab M0ller experiment is to measure the effective weak mixing 
angle, sin^ y to an absolute precision of 0.00025, or in terms of the weak charge of the 
electron, to within 2.3%. The running of sin^ 9^ in the SM gives a prediction for this 
quantity at low energies [33]. The result depends rather strongly on the scale at which jZ 
box diagrams are evaluated: Q^wW = -0.0472, while Q^(Mz) = -0.0462 [35]. Higher- 
order calculations are needed to determine the appropriate scale at which to evaluate these 
contributions. Following the suggestion of [36], we split the difference and take it as a 
contribution to the theoretical uncertainty: 

= 0.0467 ± 0.0006(Z-pole) ±0.0005(7Z) ±0.0011(experiment) (10) 
= 0.0467 ±0.0013, (11) 

where we have added the various sources of error in quadrature. Shifts from the expected 
asymmetry, similarly to Q-weak, go as 



2^2 
2 r' 

V2 
2MI,Gf 



4-4). (12) 



An interesting feature of this quantity is that it depends only on electron couplings, and not 
those of the quarks. A measurement of this quantity would break the scaling degeneracy 
that plagues the other measurements discussed previously. 



5 LHC and Q-weak Analysis 

To extract coupling information from the LHC from the measurements described in Section [3} 
we must evaluate the coefficients a, 6, and c in Eq. Q. This requires knowing the Z' mass 
from experiment. For illustrative purposes we explore the possibility M^/ = 1.5 TeV, and 
take the mass to have negiligble error [27]. Since we are not making any assumptions on how 
well the width is measured directly, we must treat this as a free parameter in our analysis, 
and fit it accordingly. This leaves five unknown parameters on which Eq. ([s]) depends: 
Ez', (iL^L-i Ql^r, urCli and dnei- Two other combinations, urCr and dneR^ appear in Eq. ([s]); 
however, these are not independent from those listed due to the following relations: 

qi^R UrCr (IrCr Cr 

= = = — = tanfc'i. (13) 

qL^L UrCl dRCL tL 
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Here we have defined the parameter 6i, which is the angle in the plane of the leptonic 
couplings. There is no way to separate q couplings from e couplings at this stage, as the Z' 
is produced via quarks and decays into leptons in this channel. This limitation is addressed 
in the next section. We include a discussion of the Q-weak measurement in this section as 
it probes the same q x e combination as the LHC 

We determine the five unknown parameters by performing a scan over this parameter 
space. For each point in the five-dimensional parameter space, we can reconstruct the pre- 
dicted LHC measurement bins discussed in Sec. |3] using Eq. ([5]), and compare to the actual 
observed values. In the absence of data, we construct measurements for the LR 



models described in Section |2.2| to test our procedure. As we present analyses for several 
coupling combinations with varying assumptions regarding integrated luminosity and energy, 
we organize our results into subsections as follows. 

• We first present a discussion of the relevant uncertainties affecting both the LHC and 
Q-weak measurements. 

• We discuss an initial three-parameter comparison between data and model hypothe- 
ses that can be performed with only a few fb^^ of data at the LHC. 

• We present results for integrated luminosities between 10 — 100 fb^^ for the Cq and 
coupling combinations defined in [8], which are most relevant when only on-peak data 
gives significant information. 

• We then discuss extraction of the q x e coupling combinations assuming large LHC 
integrated luminosities of 300 — 1000 fb~ , and demonstrate how the Z' width can be 
determined by comparing on-peak and off-peak bins. We add the expected Q-weak 
measurement to the analysis at this stage. 

• Finally, we summarize our analysis and present results showing how well all Z' param- 
eters relevant at the LHC can be determined as a function of integrated luminosity. 



5.1 LHC errors 

Our ultimate goal is to extract the underlying Z' couplings. We will need to perform a 
test between hypothesis points in the coupling parameter space and the "observed" values 
for each of our example models. For early LHC running, we will want to know whether 
one can start eliminating canonical models, and we do a pairwise comparison between our 
test models to see how well this can be done. For later LHC running, we will demonstrate 
the efficacy of our coupling extraction procedure. To do so we must study the uncertainties 
which affect our analysis procedure. We classify errors into four categories. 

• PDF . We determine the PDF errors using the CTEQ 6.6 PDF error sets [30] according 
to the prescription of [37], and account for correlations among our measurement bins. 
To be conservative, we take the larger magnitude shifts from each of the 22 eigenvector 
directions in PDF parameter space, and when forming the interpret these errors 
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as la. If the Z cross section is used as a standard candle, one should instead find 
the PDF error in the ratios of our bins to the expected number of Z events. If there 
were significant correlations between our bins and the Z cross section, this would 
reduce PDF errors. Testing indicates this is not the case, and for simplicity we do not 
normalize to Z production to determine PDF errors. 

Statistical . For each considered value of integrated luminosity we determine the count- 



ing error in each bin. The bin choices for each luminosity are described in Section 3.1 



• Theory. We have calculated all cross sections at NLO in the strong coupling. Residual 
scale errors are negligible compared to our other errors [8], and it is known that the 
NNLO corrections have a very small impact on the central value of the NLO result [9, 
38]. QED corrections in the final state may be important in predicting the invariant 
mass distribution [39,40], though with our coarse binning this issue is less important. 
Further electroweak corrections come in two forms: vertex corrections, which can be 
absorbed into the effective couplings extracted, and box corrections (where two sides of 
the box are W, Z, or Z'), which cannot. The latter have large logarithmic corrections in 
the energy, \og{s/M^^z) [41], and should be determined when examining off-peak bins. 
We assume that a full analysis will have computed these corrections; the actual shifts 
in the coefficients a, b, and c should not be large enough to affect the efficiency of our 
method. However, we do include, conservatively, a 5% error in the overall normalization 
of the off-peak bins to account for remaining errors after EW corrections have been 
included. We have found that including this error does not substantially affect the 
quahty of the coupling extraction. 

• Detector . Our binning is coarse enough that we take detector issues such as energy 
mis-measurement, cracks, etc., to be sub-dominant. We neglect them in our study. 

5.2 Q-weak errors 

The Q-weak experiment will also probe combinations of quark x lepton couplings; from 
Eq. g, it will measure the combination 

{3qLeL + 2uReL + dReL){— - 1) = V2M|,Gi.Ag^. (14) 

The Z' mass will be known from the LHC The anticipated uncertainty after combining both 
experimental and theoretical errors is large compared to the expected deviations of our test 
models: 



6 



Qw 

5AQ^ = 0.0032 



(15) 
(16) 
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For reference, we list the expected shifts due to our test models at 1.5 TeV in Table [Sj A Z' 
would have to be fairly strongly coupled to be observable by Q-weak. The Q-weak experiment 
may add somewhat to the discriminating power of the LHC, however. We assume that Q- 
weak finds a value consistent with our test model, and for each test point, we generate a 
measurement and add it to the appropriately. For simplicitly, the Z' parameters are not 
run from the Z' scale down, but compared directly. If the running is 0(10%), as expected, 
this can be safely neglected compared to the experimental error. 





AQ^ X 103 


X 


1.17 


LR 


-0.45 







X* 


2.24 



Table 5: Expected shifts in the weak charge due to a 1.5 TeV Z\ in unit of 10 ^ . This should 
be compared to the expected measurement precision of ^AQ^ = 0.0032. 



5.3 Early LHC results: comparisons 

We first assume a center-of-mass energy of only ^/s = 10 TeV, consistent with early-running 
scenarios. At low luminosities, the statistical errors are too large for a meaningful coupling 
extraction. Therefore we ask the question, is it possible to begin to discriminate models? How 
can one rule out a "model" ? High-scale models typically determine the Z' fermion charges, 
but the overall coupling running or the leptonic branching fraction may depend on other 
parameters of the theory not considered here, or even predicted by the model. Therefore, 
in this first pass to distinguish models, we wish to consider observables independent of 
these parameters. We first note the following. For very early running, the error in off- 
peak observables is too large for them to be useful. Secondly, the overall on-peak cross 
section, apeak, is proportional to the overall coupling squared and the branching fraction: 
o'peak oc g\,Br{Z' — i> e^e~). Therefore, if one takes on-peak observables normalized by the 
overall on-peak cross section, ai = ai/ apeak, they are independent of the overall coupling and 
leptonic branching fraction. Considering these observables only will allow us to differeniate 
models defined only up to these ambiguous parameters. In forming our comparison using 



the four bins described in Sec. 3.3, we normalize the bins to the on-peak cross section, 
and we restrict ourselves to on-peak bins. As these observables are no longer independent 
(^i '^i — ^peak). One bin is dropped, for three observables, and three degrees of freedom. We 
present the comparison of experimentally "found" modes versus other model hypotheses 
in Table [6j For reference, 68% confidence corresponds to = 3.5, 90% to 6.3, and 95% to 
7.8. 

In performing this test, we have assumed statistics consistent with our prior assumptions 
about the overall coupling and width for the experimentally "found" model; we remind the 
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3 fb"i 


10 fb-i 


30 fb-^ 




X ^ LR 


X LR 


X LR 


X 

LR 


- 1.7 6.7 
3.8 - 0.58 
6.8 0.3 


5.1 23.1 
11.3 - 3.3 
22.9 1.4 


15.7 54.7 
30.6 - 5.2 
55.1 2.6 



Table 6: comparison (with three degrees of freedom) between models for early LHC 
running. Hypothesis models are in columns, and are tested against the experimentally found 
row models. Errors are determined from the hypothesis model. 

reader that our procedure is independent of these assumptions. We note that the model is 
difficult to differentiate under these assumptions due to its lower production rate. We note 
that "nice" models with clean signatures such as the x model and left-right model can be 
distinguished from each other at over 90% confidence with as little as 3 fb~^. 

5.4 Early LHC results: initial parameter extractions 

For slightly higher luminosities, we can use the off-peak bins to begin to reduce the allowed 
region in coupling space. As an example, we plot the 95% confidence region for two of the 
coupling combinations (g^, x e^j and x e^) for measurements corresponding to the x model 
in Fig. [TJ If all couplings and the width are fit simultaneously, this corresponds to a 68% 
confidence region. We make the following observations. 

• At 10 fb~^, one only gets a vague notion of the size of the couplings. Since we do not 
assume the width is known from the Breit-Wigner scan at this stage, this is not due 
to knowledge of the on-peak total cross section. It is instead coming from the limited 
size of the deviations of the off-peak bins from SM predictions. 

• At 30 fb~^, we see multiple "islands" emerge; these are due to the nonlinearity of 
the cross section in the g x e parameters. The allowed coupling space has shrunk 
considerably, and we see the beginnings of a measurement for the correct island. Note, 
however, that all sign degeneracies remain in the couplings; it is still the on-peak 
measurements leading to reduction of the allowed region. To determine these signs, 
one must wait for more integrated luminosity. 

At this stage, the couphngs themselves are not measured particularly well, though we are 
seeing quite a reduction in allowed parameter space. In Ref. [8], the anticipated errors for the 
combinations of the couplings measurable on peak were derived. Given the success of that 
extraction, we wish to revisit these results using the chosen binning and lower luminosity, 
again choosing the x model as an example. We plot in Fig. |2]the 68% confidence regions for 
the four parameters Cg and (g = u,d), defined in Ref. [8] as 
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Figure 1: 95% confidence-level region for the Ql ^ gr and Ql x ei couplings extracted from 
measurements corresponding to the x model, at 10, 30, and 100 fb~^. Other model values 
for these couplings are shown for orientation. 
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2 , 2\t 2 , 2\ ^Z' 



= {qR + qi){eR + ei) 



2ATrTz' 



e. = (ql-qDiel-el)^. (17) 

These have the advantage that they are directly related to the cross section near the resonance 
peak, 

ivJlliAP " S: lanMM, + ai(MM,] (18) 

q=u,d 

where we have rescaled the coefficients a^' = "^^fy-o-i of Eq. (5) to move most of the width 
dependence into the coupling combinations Cq and Cg. For this analysis, since we use a fixed 
bin size, we also fold a factor of (1 + 2 ) from Eq. (6) into the definitions of Cq and to 
keep the results practically width-independent. This is relatively unimportant if the width 
is small relative to the on-peak bin size, or known a priori. 

The on-peak bins are doing almost all the work discriminating Cq and e^, as these quan- 
tities are directly related to the on-peak bins. The off-peak bins do not significantly help in 
reducing the allowed coupling space at these luminosities. As such, we see nearly elliptical 



errors as one would expect from a linear propagation of errors from Eq. (18). However, the 
errors differ slightly from our previous analysis [8], and we wish to mention the following 
differences. 

• These are 68% confidence regions for 2 parameters, and not Icr in each parameter. 

• ^/s = 10 TeV, not 14 TeV, which affects both statistical and PDF errors. 

• The physical constraint ul = is explicit in the parameter space scan, as are Cq > 0, 
I I *\ etc. 

• CTEQ 6.6 PDFs are used instead of CTEQ 6.5, which allows for both s quark differ- 
ences and improves high x errors [30]. 

• Less conservative assumptions are made on the correlations of differing PDF eigenvec- 
tor directions. 

We note that under these assumptions u and d coupling measurements are very anti- 
corellated on peak; the sum is known much better. This is to be expected; their separation 
depends on subtracting low and high rapidity Z' bins, thereby increasing errors, while the 
sum depends on combining bins. 
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Figure 2: 68% confidence regions for tlie coupling combinations Cq and Cg extracted from 
measurements corresponding to the x model, at 10, 30, and 100 fb"'^. Other model values 
are shown for orientation. We note that Cg and Cg also have correlations with each other that 
are not displayed. 
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5.5 Late LHC and Q-weak results 

We finally consider the improvements possible if integrated luminosities of 300 fb~^ or 1 
ab~^ are available. We now add information from a Q-weak measurement of Q^r, assuming 
it finds a value consistent with that predicted by a test model. We demonstrate in Fig. |3] 
the successive improvement upon adding different data in our extraction for 1 ab~^ at the 
LHC, by first applying the off-peak bins only to determine the couplings for the x model, 
then adding on-peak bins, and finally the Q-weak experiment. We choose to show u^eL 
versus dnei to illustrate this improvement. Due to the nonlinear nature of the cross section 
in the Z' parameters, multiple islands are possible, where one region of parameter space can 
mimic the signal of another within the allowed errors. The Q-weak experiment is helpful in 
breaking degeneracies such as these that may remain, where the false islands have a different 
enough value of AQy^^ to rule them out in conjunction with the LHC measurements. The 
difference is not in general substantial enough to rule out the islands if the LHC did not 
already disfavor them. 
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Table 7: Results of the LHC/Q-weak coupling extraction. We use the following abbreviations 
for the integrated luminosity: LL = 300 fb^^, HL = 1 ab^^. 

We use the Q-weak experiment in conjunction with the LHC at high luminosities to fit 
the couplings and width for all models. In Table [7] we list the results of the extraction of these 
parameters for integrated luminosities of 300 fb~^ and 1 ab~^. The errors listed are Icr for 
each individual parameter, while allowing the other values to float. For brevity we omit the x* 
model. Interestingly, the width is measured extremely well by the extraction. This essentially 
comes from comparing the size of on-peak bins to what would be predicted by extracting 
the couplings off-peak only; the size of the on-peak cross section is directly controlled by 
the width. We note that the precision of this measurement is not directly controlled by 
the calorimeter resolution, and is therefore controlled by different systematic errors than 
a fit to the resonance peak. The errors here are competitive with a previous detector- 
level analysis [27] of the Breit-Wigner fit with a fixed number of events corresponding to 
luminosities of 100-300 fb~^. We have tested that adding knowledge of the width from the 
Breit-Wigner fit following the study in Ref. [27] does not significantly improve the quality 
of our results. 

The couplings, especially the larger ones, are measured to good precision. The model is 
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Figure 3: Reduction in allowed coupling space with successive application of off-peak mea- 
surements only, on-peak measurements, and the Q-weak expected result for the x model, 
at 95% CL. We note the removal of the small island in the upper-right quadrant upon the 
inclusion of Q-weak. These results include 1000 fb~^ of integrated luminosity at the LHC. 
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an exception; several parameters are measured poorly. This is due to the lack of asymmetry, 
where the differences F — B are zero, giving a large relative statistical error. For non- 
leptophobic, narrow-width Z' models, we consider this a worst-case scenario. However, 
the table above does not tell the whole story. The 5D 68% CL region includes strong 
correlations between the parameters, and portrays a dramatic reduction of parameter space 
compared to a naive estimate from the table alone. We plot several 2D projections of the 
5D confidence region for our test models in Figs. |4] and [5] to demonstrate. These regions can 
also be interpreted as 95% confidence regions for the two parameters displayed in each plot 
separately, allowing the others to float. 
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Figure 4: Width versus qh^L fits for the considered models. Widths tend to correlate with 
larger couplings. 

We make a few remarks on these results. First, the x model measurement improves sub- 
stantially with increased (super-LHC) luminosity. Remaining degeneracies in the parameter 
space are removed. However, for the other models, there is only minor improvement. The 
allowed volume of coupling space is small, and our test models are easily distinguished. Also, 
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Figure 5: Coupling fits for the considered models. The actual value for each model is shown 
by the labeled black symbols. 
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we point out that it is often possible to mimic the signals of our test models within errors 
while keeping zero, as can be seen in the upper panel of Fig. |5j If the average asymmetry 
coming from the quarks is small enough, one can account for it in the leptonic couplings, 
and get the relative amount of u and d correct by adjusting the values of ur and without 
affecting other measurements. For example, the model is extremely difficult to pin down 
due to the lack of asymmetries. One cannot even favor a sign on dnCL in this model, as can 
be seen from the lower panel of Fig. [5] 



5.6 Summary: single parameter extractions 



As our analysis contained many parameters, we review here the combinations that can be 
probed by the various data sets. At low luminosity, we have seen that the majority of the 
analyzing power comes from the on-peak bin. The natural coupling combinations to utilize 



with this data set are the Cq and Cq combinations defined in Eq. (17). These quantities are 
formed from (g x e)^. At higher luminosities, the off-peak bins contain enough events to 
allow a direct measurement of the linear combinations g x e. We summarize this information 
in Table [8j Only six of seven possible parameters appear in this table; this refiects the 
g — > yg, e — > e/y degeneracy discussed earlier. 



on-peak LHC: 


Mz',Cu,Cd,eu,ed 


on-peak -|-off-peak LHC: 


Mz', Tz', QLei, QLeR, urCl, dnCi 


Q-weak: 


{3qL + 2ur + dR){en - eL)/Ml 



Table 8: Summary of which Z' coupling combinations can be probed with the various data 
sets considered. 

We summarize the analyzing power of LHC and Q-weak by displaying in Figs. [6} [7| |8| and[9] 
the improvement in the coupling extraction for the x model, one parameter at a time, as 
the luminosity is increased. For integrated luminosities of 100 fb~^ or less, we conservatively 
use a center-of-mass energy of 10 TeV; for higher luminosities we assume the full 14 TeV has 
been reached. For the Cg, Cq couplings, relative uncertainties of roughly ±50% for non-zero 
couplings are possible with 10 fb~^, while ±30% errors become possible with 30 fb^^. At 
300 fb~^ and above, the g x e combinations are directly measurable with precisions of 10% 
or better. We note that a sign degeneracy in all g x e parameters remains at integrated 
luminosities of less than 100 fb~^; this is resolved at 1 cr for all parameters at 300 fb~^. 



6 Determining Leptonic Couplings 

Our goal is to probe all Z' parameters to fully determine the underlying theory; thus far, we 
have only determined the couplings g x e. A linear collider can probe the leptonic couplings 
directly, even below the resonance peak, if the mass is known [42]. In this paper, we do not 
assume access to a linear collider, and attempt to determine how well the last parameter can 
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Figure 6: la errors for the on-peak coupling combinations Cu and q as a function of integrated 
luminosity. To the left of the dashed line, ^/s = 10 TeV, while ^/s = 14 TeV to its right. 
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Figure 7: la errors for the on-peak coupling combinations and as a function of inte- 
grated luminosity. To the left of the dashed line, a/s = 10 TeV, while a/s = 14 TeV to its 
right. 
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Figure 8: la errors for the coupling combinations qL^i^flL^R as a function of integrated 
luminosity. Note that all signs are selected at 300 fb~^, but not before. To the left of the 
dashed line, \/s = 10 TeV, while ^/s = 14 TeV to its right. 
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Figure 9: la errors for the coupling combinations UfieL,dp>eL as a function of integrated 
luminosity. Note that all signs are selected at 300 fb~^, but not before. To the left of the 
dashed line, y/s = 10 TeV, while ^/s = 14 TeV to its right. 
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be measured by other means. It is possible that observation of the Z' decaying to bottom or 
top quark final states can help break this degeneracy. We pursue here instead the alternate 
possibilities of determining the invisible width of the Z', or by using a proposed low energy 
M0ller scattering experiment. The M0ller scattering analysis in particular demonstrates a 
nice complimentarity between low and high-energy experiments. A study of the sensitivity 
of future low energy M0ller scattering experiments to Z' effects has also been considered in 
Ref. [43], although the analysis performed there did not focus on the interplay with LHC 
measurements. 



6.1 Using the width 

If we assume that the leptonic couplings are flavor-universal, we can write the Z' width in 
the limit of vanishing SM fermion masses in terms of parameters already determined, and 
the leptonic coupling: 

Tz' = ^(6gi + 'iu\ + 34 + 2ei + e^) + T^^^, 
ovr 

= ^(6(gLeL)' + 'i{uRej)i' + 3(rffiei)' + (2 + tan^ Qi)e\) + Y^ew (19) 

We have denoted the decay width of the Z' into particles other than SM fermions as r„e^, 
even if those decays are into SM modes such as W^W~ . If this extra width is determine by 
measuring additional visible and invisible decays, we can extract e| up to a 2-fold degeneracy, 
with no further input]^ 



to-V^ ± J(i^'^^r - (2 + ta,i2 emiLeL)' + 3{uneLr- + 3{dne,)') 

= 2^^, ■ '2") 

The invisible width can be probed in associated production with a Z [44] or photon [45] 
by comparing the branching fraction of invisible Z' decays to leptonic Z' decays, where the 
invisible decays due to a Z' have been identified with a missing pt cut. This would allow us 
to determine Tnew The quality of this probe is highly model-dependent, so we do not include 
it in the analysis, but point out that it could be used to determine the leptonic couplings. 

If the invisible width cannot be determined to any degree of accuracy, then in the worst 
case we have bounds for e\. For a given set of coupling combinations q x e, raising (or 
lowering) e will raise the partial width into e (or q). Since we have measured the width, 
there is a natural ceiling in either direction, corresponding to the bounds 



^i>-^^ 2^^, • (") 

, 4^fe + Jii^ik? - (2 + tari^#,)(6(,iei)= + HuneLf + iidReL?) 
2 -^-^ 2TW9^ ■ 



^If cl — > 0, one should instead solve for e/?. 
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These bounds are also highly model-dependent, as they are much looser if there is a lot of 
"extra" width. This is easily achieved in models with additional exotic particles into which 
the Z' can decay [46]. As such, we refrain from putting these bounds in our plots, though 
they are always present. We also note that if the Z' width is determined, upper limits on 
all couplings can be derived by requiring that a given partial width no more than saturates 
the measured value. 

6.2 M0ller scattering 

The future J-Lab M0ller scattering experiment may be able to observe a shift in the asym- 
metry as discussed previously. Measuring the combination eytA amounts to measuring a 
hyperbola in the — plane: cye^ = (e|. — e|)/4. Using the anticipated experimental 
and theoretical errors, 

2 2 

= 0.022TeV-^ (23) 

or, for a 1.5 TeV Z', 

S{el-el) = OmO. (24) 

The shifts due to our considered models are listed below. As with Q-weak, we ignore running 
down from the Z' scale. 







X 


-0.0436 


LR 


0.0167 







X* 


-0.0833 



Table 9: Coupling combinations contributing to M0ller scattering for our various test cases. 
None of the test models with the standard coupling induces a la deviation in M0ller. In 



Fig. 10, we plot these models with the la, 1 ab~^ LHC errors in 9i, as well as hyperbolic 
upper limits on the size of the leptonic couplings assuming no deviation from the Standard 
Model is found. The allowed regions do not actually extend to the origin; there is some 



mininum size of the leptonic couplings due to the argument of Sec. |6.1[ An overall sign in 
the couplings is unphysical, so we have chosen the sign of cl for each model for convenience 
in displaying the plot. 

Our original test models should be consistent with the Standard Model in the M0ller 
experiment for a mass of 1.5 TeV at lex, though the x model is very close to the limit. We 
therefore ask the question, if the leptonic couplings were somewhat larger, what could we see? 
As an example, we consider a test model like the previously considered x model, except with 
a larger couphng; we call this the x* model. Specifically, we replace the e/ cos 6*1^ appearing 
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Z' Leptonic Couplings, M2'=1.5 TeV 



CD 




-0.4 -0.2 0.0 0.2 0.4 



Figure 10: Measurements of the leptonic couplings from the LHC and J-Lab M0ller experi- 
ment. M0ller bounds assume no measured deviation from the Standard Model. The overall, 
unphysical, sign of the coupling is chosen for each model for convenient display of the results. 
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in models with a somewhat larger value, 1/2. Such shifts in the overall normalization 
from that predicted by the canonical breaking pattern has been considered elsewhere in the 
literature [5]. In this case, the leptonic asymmetry is large enough to produce a significant 
deviation from the standard model in the M0ller experiment. This produces upper and 



lower hyperbolic bounds in the cl - cr plane. We plot the expected results in Fig. [TT 
M0ller asymmetry gives us a hyperbola in the ei - gr plane. The LHC gives us the angle. 
Together, they pinpoint the Z' leptonic couplings to a small region in the versus en plane. 




Figure 11: Measurements of the leptonic couplings from the LHC and J-Lab M0ller exper- 
iment, for a model with larger couplings. There are upper and lower hyperbolic bounds on 
the size of the couplings. 



6.3 Putting it all together 

Once the leptonic couplings are determined, they can be input into the q x e measurement 
coming from the LHC and Q-weak to extract all Z' couplings. Explicitly, the quark charges 
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are given by 



Ql 
Ur 

dR 



UrCl 
dRCL 



(25) 



We list in Table [TO] the final expected results after solving these equations for the quark 
charges and propagating through the anticipated M0ller errors. We include all LHC data at 
1 ab~^, and the Q-weak and M0ller experiments. Each parameter range given is 68% CL. 
Note that for typical models, we only get upper bounds on the sizes of the leptonic couplings, 
and lower bounds for the sizes of the quark couplings. The signs are determined from the 
LHC. We note that the bounds derivable from width constraints discussed above have not 
been included in this table. 





X 


LR 




X* 




to 0.323 


to 0.228 


+ 


0.306±°:?^^ 


eR 


to 0.111 


-0.343 to 




o.io2±°:°i 


Ql 


< -0.051 


< -0.009 


+ 


-0 102±°°20 


Ur 


> 0.045 


> 0.126 




0.102±0;057 


dR 


< -0.150 


< 0.162 


? 


-0.306±0:J^? 



Table 10: Final determinations of the couplings using data from the LHC, Q-weak, and low 
energy M0ller scattering. The results assume 1 ab~^ of integrated luminosity at the LHC. 
The sign of is chosen to be positive to fix the overall sign degeneracy, rf/j remains totally 
unknown for the ijj model extraction, and for other parameters only the sign is determined. 



7 Conclusions 

We have presented a thorough look at the ability of the LHC to analyze the underlying model 
parameters of a new Z' boson in the dilepton channel. In addition, we have studied the 
complimentary information that can be obtained from upcoming high-precision, low-energy 
experiments. Our analysis is as model-independent as possible, and includes the major 
sources of error which hinder the parameter extraction. We have discussed the important 
role played by off-peak data and the future Q-weak measurement of the proton weak charge 
in removing sign degeneracies that remain if only LHC data on the Z' peak is analyzed. 
We also have shown how a precision measurement of low-energy M0ller scattering can break 
the q yq,e e/y scaling degeneracy inherent to the dilepton channel at the LHC, and 
permit an extraction of the individual Z' charges. Our study demonstrates that a global 



32 



analysis of data from many different observations is a necessary component of any attempt 
to determine Z' couplings from the bottom up. 

Our study has elucidated the how precisely Z' couplings can be measured as a function of 
integrated luminosity. Our LHC analysis relies upon appropriate binning of dilepton events 
in rapidity and scattering angle: u and d couplings are correlated with higher and lower Z' 
rapidities respectively, while left- and right-handed couplings yield different lepton angular 
distributions. As the parameters extracted and the bins used vary as a function of integrated 
LHC lumunosity, we summarize below how we believe the analysis of a dilepton signal at 
the LHC should proceed. This template assumes that the resonance has a mass, width and 
total cross section roughly in the ranges considered here. For significantly different values, 
the timeline presented below must be adjusted. 

• With less than 10 fb~^, the first property to establish is the spin of the observed 
resonance [19,20]. At this point, a comparison between data and model hypotheses 



based on the three on-peak observables discussed in Section |5.3| allows initial guesses 
as to the underlying Z' model to be tested. This analysis is less sensitive to additional 
decay modes and overall coupling than comparisons of the total cross sections. 

With 10 — 30 fb~^, enough on-peak data is available to attempt an extraction of 
the couplings. Sufficient data exists only on-peak at this point, indicating that the 



natural parameters to study are the Cg, defined in Eq. (17). This analysis can be 
performed in a simple, model-independent fashion by following the procedure described 
in Section 5.4 We have found that the Cq coupling combinations can be measured with 



relative uncertainties of ±30% for the Mz' = 1.5 TeV test cases considered. 

On-peak observables are sensitive to only the combinations (g x e)^. The sign ambigu- 
ities can be removed by studying off-peak data, which probes the q x e combinations 
directly. This becomes possible with 100 fb~^ of integrated luminosity. We have found 
that qxe combinations can be determined with relative uncertainties of approximately 
10 to 100%, depending on the particular model and coupling combination, for luminosi- 
ties of 300 — 1000 fb^^. The couplings display very strong correlations, which indicate a 
much stronger reduction in parameter space than the individual errors would suggest. 
A precision measurement from Q-weak helps remove sign degeneracies that remain 
after LHC running. Another interesting aspect we have seen in our study is that the 
width can be determined precisely at this stage without a direct measurement of the 
Breit-Wigner shape. 

Some direct measurement of either lepton or quark couplings is crucial to break the 
q — >• yq,e e/y degeneracy that exists in the LHC dilepton mode. We have found 
that a proposed M0ller experiment can determine the combination e\ — e| directly 
in some situations, and bound it in others. Measuring this last parameter can tell us 
all the Z' couplings to Standard Model fermions, by separating out the well-measured 
combinations qxe from the LHC. 
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Our study can be extended in several interesting ways. It would be interesting to apply 
our analysis procedure to other Z' models to more fully determine its robustness. Another 
possible extension would be to add other LHC final states such as tt, which could also 
potentially break the scaling degeneracy from which the dilepton channel suffers. We look 
forward to studying these possibihties and hope that Nature provides on opportunity to 
apply our procedure to LHC data. 
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